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ABSTRACT 


A brief survey of sone of tlie methods of improving angular 
resolution of a VHP ground based radar is carried out. The 
concept of spatial frequency is highlighted, A scheme analogous 
to holography is proposed using the idea of a DPT on spatial 
samples to obtain a spatial frequency spectrum. A DPT implemen- 
tation using table look up instead of a hardwaro multiplier is 
suggested and simiiLated, Elementary smaulation of the scheme is 


carried out 



CHAPTER 1 


INTRODUOTION 

A grou3ad~based surveillance radar used, for tiie detection 
of aircraft has to, in addition to detecting the presence of 
aircraft, be able to determine the position in terms of slant 
i-ange, azimuth and elevation angles. Such position has to be 
determined to the greatest possible accuracy. In addition to 
locating the position accurately it sho-uld also be possible to 
differentiate between two targets proiimate to each other. The 
closer the proximity the better the resolution. Betcer accuracy 
does not necessarily mean better resolution nor vice-versa. 

In this thesis we are concerned only with obtaining 
improved angular resolution in the aamutli plane. If corres- 
ponding resolution is not available in the elevation dimension 
also^ the system will not be useful for posi'cioning friendly 
aircraft in an advantageous position vis-a-vis hostile aircraft. 
It could of course be helpful in determining the number of 
hostile aircraft, unless they were ’ stached' one above the other. 
As a navigation aid to friendly aircfaft igs usefulness would of 
course be enhanced, since aircraft normally know their attitude, 
their range and azimuth angle from a reference point can be 
passed to them with a better accuracy. 
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The noriaal means of achieving this has been to use 
tiigher radiated frequencies so that for given mechanical load 
(for aerial rotation) and structural complexity narrower beams 
and hence better resolution is achieved. In this thesis we 
consider a system using a radiated frequency in the VHP range. 

Use of VHP frequency confers two advantages : 

a) freedom from weather clutter and'angels' [1]. At 200 MHz 
the returns from rain is 60 dB lower as compared to that 
at 10 GHz. 

b) for a given pulse repitition frequency (p.r.f) a lower 
transmitted frequency (larger wavelength) gives a higher 
blind speed. As will be seen later for the case under 
consideration, the blind speed is far above the normal 
operating speeds of aircraft. 

The thesis is organized in the following manner. Chapter 2 
deals with some methods which have been px’oposed/used for impro- 
ving angular resolution. The chapter after dealing with optical 
holography attempts to draw a picture of an analogous system in 
the microwave region. This leads, m Chapter 3, to specifications 
of the solution to the problem at hand, based on the holographic 
analogy. The chapter spells out the details of the solution. 
Chapter 4 discusses some of the routines used in the simulation 
as well as consideration of SNR and clutter to signal ratios. 
Chapter 5 gives the conclusions arrived at as well as some 
suggestions . 
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CHilPTBR 2 

METHODS POR IMPROVING ANOULA.R RESOLUTION 

Before proceeding to details of the proposed scheme, 
a brief review of some of the methods used for obtaining 
improved angular resolution is presented below. 

The methods considered are : 

a) Superdirective arrays 

b) Interferometric methods 

c) Frequency swept imaging 

d) Synthetic aperture methods 

e) Microwave holography 

2.1 SUPERDIRECTIVE ARRIYS 

A superdirective antenna is one having a beamwidth (in 
radians) c..»c» l/largest dimension of the array in wavelengths 
[2], A characteristic feature being their dependence 
on heavy cancellation between closely spaced radiators. 

Thus requiring a high* precision in feed and placement of the 
elements. Thus not only the achievable gain of such arrays 
IS to be considerod. in addition one must keep in view the 
sensitivity of the lobe pattern to errors in feed and 
placement. As mentioned in [2] these arrays are high Q small 
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bandvridth. devices^as the amounc of reactive energy in the near 
field IS very large compared to the energy radiated per cycle, 
Schelkunoff [3] has also shown that though directive gam can 
he incx-eased the efficiency and bandwidth reduces. 

hewman et al [4] q.uote results obtained for a 9 element 
Imear broadside array of copper halfwave dipoles with overall 
array length of /4 for ob bainmg a directivity 8.5 times ’chat 
of a single dipole, fhe efficiency was 10~^^ percent and 
accuracy required was one part in 10^^. fney have however, shown 
that for a receiver array where the system performance is limited 
by external noise the poor efficiency does not effect the 
signal to noise ratio (provided (antenna efficiency x external 
noise)3?‘> internal + receiver noise). Such a system cannot obvi- 
ously be used in the situation under consideration, 

2.2 HJlBEPBROiyiEfRIC D-IETHODS 

By the term ' mterf erometric’ methods we mean those 
methods where the antenna aelements are spaced apart far enough 
to give grating lobes, i,e, if the output from the elements 
were additively combined. ¥e consider two methods. 

2.2.1 1st Method 

This method was proposed by M, I, Skolnik [5], The method 
IS based on recogniging the change in grating lobe pattern 
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that results from a change xn^l /d rabio (d is element spacing). 

He considers a radar covering a volume of space with one "broad 

transmitting beam and many fixed narrow receiving beams. The 

receiving array and transmitting antenna are separate. The 

receiving array consists of 11 isotropic elements spaced 5 ,75^1; 

apart, giving a main lobe of 0.9*^. The grating lobes are at 

+ 10°, + 20,4°, + 31.5°, ± 44.1°, + 60,5°. The region between 

the grating lobes is covered by 9 additional beams genera,ted 

by beam forming networks, Ey transmitting on another freiuency 

(or alternately using another array with different d) the 

position of the grating lobes can be made to shift while the 

main beam remains in the same position. By talcing the new 

d/A ratio as 5,2 a shift of one beamwidth in the 1st grating 

lobe IS produced (the 2nd grating lobe will shift by twice, 

the 3i*d by three times and so on). Thus ambiguities can be 

resolved. This method appears to be suited for the problem 

under consideration. Using two sots of antenna with different 

inter-eloment spacing as above one could expect a feasible 

solution. In the presence of a number of targets and clutter thi'"' 
might be difficult, 

2.2,2 2nd Method 

This method is the multipT,icative or correlation array, 

A two-element correlation array is shown in Big. 1, 
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The output is a function of the relative delays (between the 
output of the antenna element) and hence of the angle of 
arrival. The output given by 


where 


D^^iCt) = I Y(t) Y(t+T;) -dt 


X = 2nd sin 0/ C 


2T 1 Time period of received sinusoidal waveform. 


The above expression closely resembles the autocorrelation function 


‘ = lim ^ . V(t) Y(t+T) dt 

In most cases is proportional to /'pp(''^)* 

Dii(t:) represents the directional response of the array (as 
v is a function of 0). Since is also proportional to the 

autocorrelation of the received waveform it gives the fourier 
transform of power spectral density of the received waveform. 
Thus, if we use a waveform with a autocorrelation equal to 
6(0) (e.g, white noise) the angular response is similar to a 
delta function at angle Q = 0, If we used a pure sinusoid the 
directional response would exhibit grating lobes. 
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The above is a manifestation of tne space frequency 
equivalence [6]. This equivalence concept indicates that 
two widely spaced receivers of small size can achieve a 
directivity against a multifrequency source equal to the 
directivity obtained from a fully extended linear array of re- 
ceivers of twice the length against a single frequency source. 

It must be noted that no mention is made of the received 
signal to noise ratio, 

Pigure 2 taken from EOCKS paper illustrates this 
equivalence. 

McCartney [7] concludes that as with many processing 
systems employing nonlinear devices the signal/noise performance 
IS poorer at low SNR. Thus in radar it would mean reduced 
range. However, the range resolution is also poor. This is 
because pulse widths have to be large enough so that for the 
lowest modulating frequency two relatively delayed signals 
will always have an integral number of cycles of overlap, 

Further though the directional response is narrow the presence 
of other targets affects it thus the resolution is poor unless 
the returns from the second obQect is uncorrelated to the first 
objects response - a condition which cannot be expected in a 


radar 



VARIOUS SPACE COMPLEXITIES A l-ELEMEHT CORRELATI{»« 

OF A RECEIVING ARRAY SIMILAR BEAM RECEIVER OPERATING AGAINST 

OPERATING AGAINST A SINGLE PATTERNS FOR A SOURCE ^-AVNu VARIOUS 

FREQUENCY SOURCE BOTH CASES F=ECUShCY COLPLSX!T*£S 
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Fig, 27 — ^An equivalence exists in a receiving array between the coiuplex:it> of its spi'ce conhguration and the 

complexitv of the frequency configuratior <>! the source 

Fl(^ - £ 
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2.3 lEEQUBNCY SWEPT imGIKG 

Proposed by Earliat [8]. As stated by him “ 'An essential 
part of any coherent imaging system is an imaging aperture 
that furnishes a spatial sample of the wavefield scattered by a 
coherently illuminated oboect'. Images are then formed from 
the sampled data thro'ugh basically a fourior or Presnel transform 
operation. 

Two coherent xiravofronts will interfere and cause interference 
patterns just as in the optical case. The concept of frequency 
Swept imaging is based on the 'breathing' (expansion or contra- 
ction) of the pattern when the illuminating frequency is changed. 

Consider Fig, 3. 01 and 02 are scattered plane waves of 
unit amplitude, R is the reference plane wave with amplitude 
A 1, The time varying component of tne receiver output at 
B IS 

I(t) = 2A[cos k^B(sin - sin ^qjj) (1 + 

+ cos -j k 2 _B(sin ©Q 2 - sin + at)( ] (2,1) 

where a linear frequency sweep w(t) = w^(l+at) of duration T 
IS used where w^ is the initial angular frequency. This may also 
be written as 


A = ''■pCl+at) 


( 2 . 2 ) 
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If /, 2 is tile terminal wavelength., then, 

aT = (/I J_/,, 2 - 1) (2.2) 

The above expression can be put as 

jan(l+at) -oa. ( 1 +at) e 

I(t) = 2A[(® ± 2 ^-- — ) + ■ 

where 2 (sm Qq^ 2 ■*■ 

3 a, (l+at) 

Considering the fourier transform of one term e 

3 a., 3 a, at - 3 wt 

P(w) = I [e • e ] * s • dt 

3a^ 

= e • 6 (w - a^a) 

However, since t can only take values such that 0 t T, 
The above transform has to be convolved with [Tsinc(wT/2) ], 

¥e, therefore, have F(w)= e T sincL(w - a^^a) T/2] 

Similar expressions can be obtained for the other terms. 
Therefore the fourier transform of I(t) is 

3 a, - 3 a., 

I(w)s=ATLe sine ^ [w-a^^a) T/2 + e sine- (w+a^a)T/2 

+ e sine (¥-'a 2 a)T/ 2 ' + e sine (w+a 2 a)T /2 ] ( 2 , 4 ) 
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This now represents the image reoons'cruction along with 
bhe conjugate image parts. By the Bayleigh resolution criterion 
the minimum resolution is obtained for 


aitt - a^a = ^ 


Therefore, a(a^ -- 3 ^ 2 ^~ 


A-n 

T 


a(k^ sin sin 


Substituting for a 


411 

T 



(sin ©Q^ - sin ©q2^ 



Substituting 

sin - sin 0^2 = [ ; ; - 7-- ^4 t] • s (2.5) 

If we had a sweep of 2 ~ 18 GHz and a B = 1 km the resolution 
for small values of ©q^^ and 0 q 2 as worked out below. 

^ = 150 X lO"^ m /■ 2 = ^ 



—6 

150 X 10~^2. radian 
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Tills IS a resolution of metres at 100 kms. 

This method again is not suitable because it involves 
wideband transmitting systems, 

2.4 SYNTHETIC APERTURE 

In this method a large aportui-e is synthesized by 
sequentially transmitting (and receiving) from different 
points. This in effect amounts to time sharing a single 
transmitter/receiver by the elements of an array. The returns 
from each point are stored] weighted and added. The principle 
has been used for airborne synthetic aperture radars used for 
terrain mapping. Synthetic apertures are also used in radio 
astronomy. In this the motion of the earth provides the 
necessary translation of the receiver. An example of such 
a system is the Ooty Radio Telescope worhing at 327 MLz [9]. 

Harmuth [10,11] discusses stationary synthetic aperture 
radars based on nonsinusoidal functions. 

The similarity between the coherent processing of the 
airborne SAR and holography has often been highlighted [12,13,14,15]. 
As pointed out by Koch [14] a feature which sets holography and 
coherent sideloohing radar apart from ordinary photography and 
radar is that it is not 3 ust one plane that is in focus. This 
IS because each point forms its own zone plate. Since formation 
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of holograms do not require a doppler shift such doppler free 
coherent radars should also be possible. ‘ They could take 
tne form of long linear arrays of independent receivers, 
crossed arrays (Mills Gross) or square arrays’ , 

¥e now burn to Microwave Holography, 

2.5 MIGH0¥AVB HOLOGRAPHT 

In conventional photography the photographic film records 
the light intensity over the exposed portion. The incident 
light however has an associated, phase. The variation of phase 
(which depends on the wavelength and the angle of arrival) over 
the film however is not recorded. This variation of phase 
over the filia if the illumiiiation is coherent alongwith the 
mtensity variation provide the total information of the 
incident light wave. It may therefore be used to reconstitute 
the incident wavefront. 

By introducing another coherent wavefront called a 
reference wave, we obtain an interference pattern which reflects 
the phase information. Consider fig, 4. 

In Big, 4(a) we have only illumination by the coherent 
incident wave of wavelength /\ of amplitude A^and a phase 
variation along x of ^ x sin 0, (We consider only the one 
dimensional case which can readily be extended to two dimensions). 
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J'lg. 4 


P’a -t ©graphic film 
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The film will therefore have a uniform exposure proportional to 

2 

the intensity A. 

In. Pig, 4(h) a conerent reference wave has been introduced. 
The reference wave since it is normal to the film has a constant 
phase along The reference wave has unit amplitude. 

The two wavefronts will nox^ produce maxima and minima 
depending upon xfhere they are in phase and where in antiphase. 

In the above case bhe ma^imas will be seperated by 

X = /sin 0 

Thus the variation of intensity over the film records the phase 
information. The amplitude of iLicident xrave determines the 
amplitude of the variation i,e, bhe iiitensity of the maximas and 
minimas. Since only the real part will bo recorded the film 
will record a pattern given by 

A oos (^s.nyi-2) fzo 

1 A 

Vith reference to Pig, 4(a) we have mentioned that the 

phase (cp) as a function of x is given by (considering only the 
X dependent part of the intensity), 

) X sin © (taking (p = 0 at x = 0) 

i,e. d(p = ) sin © dx 

■ k il = 


(2.6a) 
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The left hand side is termed the spatial freluenoy 
(f^). ViHiere the temporal frequency has units of cycles per 
unit time, spatial frequency has units of cycles for unit 
distance, 

1/e see that 

f^ = sin 0/^ = l/distance "between interference fringes 

Now if we took - 0 as the angle of incidence the spatial 
frequency would be - sin 0/_^ but the spacing between the fringes wo 
would still be the same, ¥e thus have an aaibiguity an the phase 
information 'coded’ in the fringe pattern. It is easy to see that 
this arises from the fact that only iiitensities are recorded. 

Thus the phasor which leads the reference phasor by a waJLl produce 
the same resultant amplitude as a phasor which lags by a, 

a Now the exposed film is illuminated by a coherent 

ctrf\j>Li t-ucies 

plane wave and the lat ens i - t4: -&e observed in a plane parallel to 
the plane of the film and situated at a distance Z from it as 
shown in the Eig. 5. 

An elemental area in the xy plane will transmit amplitude 

■^jl_(2!:,y) dx’dy. I'/hen observed from point; Z in the ZI plane 

this results in an ainplitud,e l_A^(x,y) dX’dy]/Z (assuming 

Z^ y (x-rZ)® + (y+T)*), The associated phase change is 

) V(Z+x)® + (T+y)^ -f Z^, The resultant at Z of all such 

elemental areas is given (the discussion does not consider the 
part of the exposure independent of (x,y)i. The amplitude 
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transmitt^^^ is proportional to the exposure [a). 
A 


f I ^ • C ^ f 5^ ) 

iX,J) = ; j — exp[o(2it/^ ) V'U^zJ^^TiStP+Z® ] • dx dy (2.7) 


Sxmplifying 


iTCZ + x )* +“( Y + y )*“'+ 2 * = ZU +(^)2 + ( 2 ^)* 


Taking the first terms of the binomial expansion 




Tresnel Approximation 


= Z + [(X+x)* + (Y+y)*] 

' h,(x,y) _ 2^ 


Ao(X,T) = i 


il 


2 exp[-3 ‘Z] expj_--o ^ ^ 


[ (Z+x)* + (Y+y)®]] dx dy (2*0) 

exp[-o ^ Z]- exp[-3 X Iz^X^+Y^)] L cti i / a a>T 

a jj A^(x,y)exp [-3 +y )] 

• exp [-3 (xX+yY)]- dx dy (2.9) 


making the Fraunhofers approximation Z r- x^ + y* 


J exp,-, . y(^)]3 ax ay 


¥e put fj = X/(Z ), f^ = Y/(Z ) 


(2.10) 


Since the film had recorded the intensity as given in eqn, 2,6 
the intensity of the light transmitted mil be given by 
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A^(x,y) = cos ( 2 % + B (2ai) 

■where B is the constant ■backgro^md illumination. 

Therefore, the term under the integral sign now "becomes, 



+ B 


cos (21!; f 

/ 


X 

exp [-3 


x) exp [-3 

271 X 


27t fy x] dx • I exp [“3 

.A. ^ 

-3 2it y] dy dx 


271 fy y] dy 

( 2 . 12 ) 


If we took the limits of the integration as +'' ' we see 
that both terms are fourier transforms. "Using the fourier 
transform relation for a sinusoid and a constant we get the 
result as 


- ^x^ ^x^3 * + B6{fy) . 6(fy) (2*13) 

Thus we would have got three bright spots on the X-axis, One at the 
the origin due to the background illumination , The other two at 

= ± ~ ®)/?l (2.14) 

or X = + Z Sin 

However, because of the Praunhofer approximation we cannot 
take the limits of integration as (also because of the 
finite size of the film). To account for this we can as usual 
multiply A^(x,y) with a rectangular window function correspond- 
ing to the X— dimension of the film (as we are considering only 
the one dimensional case). Thus the delta f'unctions are 
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convolved with a sine fimction. Because of this, along with 
the bright spots at Z = 0, + Z sin 0, we would get secondaiy 
bright spots caused by the sme function side lobes. 

Above we have seen a non-dim ensional case of the concept 
of spatial freq.uency and how the recording of the fringe pattern 
and subsequent reconstruction are intimately related to it as 
also the cause of ambiguity in reconstruction. 


However, in general the iOonination will be complex with 
amplitude and phase variation over the surface. Since the film 
recording will depend only upon the real part, i.e, 

A. (x,y) + A?(x,y) 


Substituting this for A. cos — ) in eqn. (2,12), 

A 


¥e have fir the first term 


h’Ai(x,y) + A*(x,y)] exp[- 3 2 'n;f 2 -x]dx | exp[-o Zuf^y ]dy (2,15) 


If the reference illumination is changed to an angle a the 


reference itself would have a phase variation which would 
produce an equal change in the phase of A^(x,y) which would 
now be 


Aj^Cxjy) exp(o2'iif^x) 


The above integral now become 



( 2 ^ 16 ) 


J 


exp[-o 27jfy.y jdy 


I 

‘ j [Aj_<x) exp[-3 27tx(f^-f^) 


+ A*(x,y) exp[- 32 Ttx[f^+f^]] dx 


Here if A.(x,y) were caused Toy diverging beam-A?{x,y) would 
represent a converging wavefront. Further the spatial spectrum 
on the Zo»axis iS now seperated. The virtual image is about 
X = +f^ and the real image about X = -f^* 

u ^ OC 


The above is a one~dimensional form of the well-known 
Leith-Upatniek hologram. 


It must however be remembered that the fourier transform 
relation between the amplitude distribution along X and the 
film exposure along x is valid only if the Fraunhofer condi- 
tion IS satisfied. We have also seen how under this condition 
the optical reconstruction is equivalent bo a fourier transfor- 
mation, The above is a simple one dimensional formulation of the 
well known result, that in the Fraunhofer zone the field distri- 
bution IS a fourier transform of the aperture distribution [16], 

The fourier transformation is between space domain (x,y) and 
spatial frequency domain (f^,fy). The above discussion is as 
applicable at radio and acoustic frequencies as at optical. We 
have seen that the spatial frequency represents the angle of 
arrival of a wavefront at the plane of aobservation. Therefore, 
their spectrum will lie between + l/^ . Therefore, the maximum 
required spatial sampling rate is 4/^] samples/unit distance. 
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However one practical dilficulty at radio frequencies 
compared to optical frequencies is th.e large size of the 
holographic 'plate' further corresponding to a photographic 
plate no good area detectors exist for radio frequencies. 
Experiments done at 35 MHz using liquid crystal area detectors 
required power levels at the detector of the order of a few 
milliwatts/ cm* Ll'7]- Such high signal level would he considered 
a radiation hazard'. Aircraft returns cannot he expected to he 
of this order. 
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GHAP?flR 3 

STATEmm OP THE PROBLEM AHD SUGGESTED 

SOLUTION 


The problem taken up for study is the xmprovement of angular 
resolution for a groimd based radar operating in the vhf range 
(100 ~ 300 Mc/s) and used for the detection of aircraft. 

Some of the representa bive parameters of the radar set 


are given below, 

a Peak transmitted power 

b Puls e width 

c Transmitted frequency 

d Ifex detection range on 

(on 1 m^ target) 

e Pulse repitition frequency 

f Antenna scan rate 

g Transmitted beamwidth 


750 m 

6 [IS 

150 MHz 
350 Kins 

400 pps (giving an unambiguous 
range of 375 1ms) 

6 rpm 

15° in azimuth 60° in 
elevation 


Further antenna size increase for improvement of resolution 
IS not considered feasible. 


The Scheme Envisaged ; 

The transmitter system of the available radar is retained. 
The receiver system however would consist of one or more antenna 
anrays and receiver processors situated at a distance of 



26 


500 lOOO metres from the transmitter site. Some of the 
operational advantages that could result from such a scheme 
are ; 

a) The transmitter and its antenna system need only "be 
a cheap 'hare-bones' affair thus easily replaced if 
destroyed. Thus reducing system vulnerability to anti- 
radiation missiles. 

b) By giving operational cover at lower frequencies also 
it will force a potential jammer to expand his energy 
over a wider spectrum. 

c) AS the transmitter and its antenna could be relatively 
simpler and cheaper deployment of decoys would be much 
easier* 

The price to be paid for all this will be increased pro- 
cessing complexity and a sophisticated data cOiimunication system 
to link tne receiver sites. 

Approach to a Solution 

Prom what we have seen iii Chapter 2, we now make the 
following observation, 

a) Our problem is then to determine the spatial frequency 
spectrum as this will give us the different directions 


of arrival 
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b) The spatial spectral bandwidth can have a maximum possible 

value of 2/^ , i.e.[(cos to cos {%l2)/^\ ] 

c) By invoking the sampling theorem it should be possible 
to sample only at discrete spatial points and be able to 
reconstruct the original spatial information. 

d) As is well known a convenient way to do this would be 
by doing a DFT using the FPT algorithm-. 

A number of experiments have been done using the holo- 
graphic techniques and then sampling the interference pattern 

The method of forming interference patterns to determine 
the phase though invaluable at optical frequencies is not 
quite appropriate at radio and acoustic frequencies as well 
established methods for coherent detection already exist [18]. 

We need not, therefore, have recourse to physically generating 
interference. However, the notion of discrete spatial samples 
of both the amplitude and phase of the spatial information 
and the notion of a spectral analysis /filtering of this discrete 
sequence could possibly be used. 

We are forming a beam by processing. By apx)ropriate 
filtering the 'beam* may be 'pointed' in any direction at any 
time after the spatial samples have been recorded in much the 
same manner as after a holographic record has been made we may 
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view the scene by illuminating the hologram*, the view depending 
upon the position occupied by che viewer. In our case now the 
DFl corresponds to the illumination and the position of the 
viewer corresponds to the frequency we are interested in. 

At this point, it would not be out of place to mention 
a typo of frequency scanning technique known as within the 
pulse scanning. In this system the relative phase shifts 
of the array elements are varied by rapidly varying the phase 
shifts of the LO frequency fed to each receiver element [19] . 

One of the systems reported [20] for scanning is given in 
Fig, 6 , But, as no record is kept of the waveform received 

at each receiver, on phase shifting and combining, information 
from directions other than where the 'processed' beam points is 
lost. Because of bhis sampling caused by within the pulse 
scanning only a fraction of the energy returned fx''om a target is 
utilized. This loss of signal energy will be a disadvantage. 

In the proposed system however the scanning is not done 
at EF or IT but after coherent detection. 

To summarize the above, : It is proposed that the system 
use a array, the outputs of all elements of the array be 
coherently demodulated to obtain a complex quantity (the In-phase 
and Quadrature voltages). Analogous to the holographic method 
a Dl'T be done on this sequence of complex samples to obtain 
the discrete spabial frequency spectrum. 
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We now tiirn to the design of such a system. 

3.5 Do tails 

We have seen above that ohe spatial freq.uency bandwidth 
being 2/^ the minimum sampling rate would have to be 4//\ 
samples/unit distance. It is here that we take advantage of the 
bandlimiting done by the t2?ansmitter beam. Because the trans- 
mitter beam does not extend from -ti/B to 'n;/2, but has a beamwidth 
of 15 '^, as mentioned, we will have a spatial spectrum correspond- 
ing to this 15 ° only. 

As the rate of change of sin 0 (w.r.t. 0) is maximiua when 
© = 0 the maximum spatial bandwidth will bo obtained when the 
beam lies from (-7*5° to + 7*5°), The spatial bandwidth now 
will be 0.26/y^ requiring a sampling rate of 0,52/^j samples/ 
unit distance* 

The assumptions in the above are : 

a) The transmitter sidelobes are neglected 

b) The targets of interest lie at ranges greater than 30 Ems, 
Thus with a separation of 0,5 ICms between fcne tiansmitter 
and receiver systems the angular co-brdinate to a targev 
from either point is same. 

If we remove assxnnption (b) all it would mean would be 
that the locus of equal delay, instead of being a circle, is 
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an. ellipse with the two sites as the focii. 

We have now obtained the reqtuirement that the antenna 
should be seperated by a distance of /j /0,52 (sampling rate 
0.52/^ sample/imit distance). 

Now we would like to know how many of such antennae are 
required. Ihe more number of samples that we obtain, the better 
the resolution. As the number of such samples is finite, 
consider it as applying a rectangular window to an infinite 
sequence. The infinite spatial sequence if transformed would 
give a delta function in the spatial frequency domain (this is 
the ideal case). The rectangular window however, convolves a 
sine function with this. The transform will now be a sine 
function given by 

n sin( 27if Y D/2 ) 

where the spatial samples have been taken over + D/2, f^ is the 
spatial frequency. The first nulls will be at 

27tf^ D/2 =+71 [or - f ) D/2 = + ^ if we consider 

resolutl:on about some other frequency f^] 

f ^ = ± 1/D 
sin 6' = + / /D 

= + sin“^(/) /D) 
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Thus for two ohoects to be resolved we use the Hayleigh 
criterion that they must be separated by at least 20’ . 

If we took 0* = 1°, i.e,, two targets separated by 2° 
should be resolved then, 

D = 57.2/ 

The number of spatial samples = 57. 2^^ x / 30 samples 

Since we propose to do a PIT on these samples we take the 
number of samples as 32, Thus the number of antenna used is 32, 

How for 32 samples we determine what would be the resolution 
at an angle of 0° and at say 60° (since the resolution will get 
poorer as we approach the endfire position of the array). 

at broadside 

O' = sin"^[/' /(32x-i /0.52)] = siir^L0*0162] 

= 0,95° 

Therefore, resolution = 1.86"' 

At 60° 

[sm 0^ - sin 60°] = + ^ /(32x ,! /0,52} = + 0,0162 

Sin ©2, “ i 0,0162 + 0,865 

©2 = (0.8488), sin“^(0.8812) 

4, 58.3°, 62° 
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The beamwidth/^is approximately 1«8° which gives a 
resolution of approximately 3.6°, 

biiiilarly, at 45° the b earn width ^ is approximately 1.3*^ giving 
a resolution of 2.6°, 

Thus we See that even at 60° we have obtained a satis- 
factory reduction of the transmitted beamwidth of 15°, 

In a normal surveillance radar two targets at the same 
range and azimuth angle but different heights cannot be resolved. 
In the system that we are considering, targets which lie at the 
same range and lying on the surface of a cone with the linear 

Uiiil r\tb te 

array as its axis^ In other words, wnereas in a convGiiti.onal 

system the co-ordinate system consists of the length of the 

and 

position vector (range), /the angle of the position vector 

projected onto the horizontal plane^ from a vertical plane 

passing through the reference direction (North) this is the azimut 

angle - in the present system the angle of the position vector 

axis 

from the linear array /is whr.t is obtained. This is shown in 
Pig. 7. 

As the transmitter beam rotates in azimuth from North (000) 
clockwise to Bast (090), South (180) and West (270) what is the 
corresponding behaviour of the angle output at the receiver 
system after processing*? 
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The DFT will give us 52 spatial frequencies over a 
spectrum of 0 to 0,52/^^ , Are the different sectors correspond- 
ing to this spectrum distinct or do they overlap? 

The sectors (0 to 30°), (30° to 90°), (-30° to 0°) and 
(-90° to -30°) will overlap (by approximating 0,52/y to 0,5/^> 
otherwise instead of 30° it would be 31.4°). These angles are with 
respect to the array axis. Therefore, they correspond to 8 
sectors in the azimuth plane. One way that this ambiguity can 
be resolved is by having the transmitter antenna rotational 
information available also. 

Thovhf antennas will have to have RF amplifiers mounted on 

them to boost the signal level. Such an antenna amplifier 

has been described in [21], After mixing with a stable local 

stages , 

oscillator and D? amplif ication^whose gains are adjusted to 
obbain necessary taper, coherent detection is performed by 
multiplying the signal with the inphase and quadrature outputs 
of an oscillator phase locked to the transmitter pulse. After 
lowpass filtering we obtain the inphase and quadrature compo- 
nents of the signal received. 

We have now obtained 32 spatial samples. Each sample 
consisting of a real and imaginary part. These must now the 
be digitized. Since the bransmitter pulse width is 6 ps bhe 
range gate will also be of 6 \is duration. Thus bhese 32 samples 
must be processed every 6 ps. If we want a slack of 1 ps then 
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the conversxon must be done in 5 hS, The AD 574 A/D converter 

approx. 

has a conversion time of/So \is for 8 bits. Thus 6 A/D convertors 
each for real and imaginary parts would be required at each of 
the 52 receiver outputs giving a tobal requirement of 384 A/D 
convertors. These 6 A/D convertors would each deal with every 
sixth range cell, Tne receiver output being connected to each 
in a round-robin fashion^ The outpu'c of the A-D convertor would 
have to be switched similarly to maintain the range cell sequence. 

We now have 32 samples each of 16 bits (real and imaginary) 
these will be updated every 5 ps. A DIT is performed on these 
samples. The algorithm used is the one given in Pig. •§ ^ [22]. 

Base in coding a simulator for sa,mples which are a power of 2 
in number was the reason for the choice, Por each of the 5 
stages we could have 32x16^ bit locations for storing the 
results of each stage. In which case our requirement is that 
each stage be finished in 5 [is , In cose wo use 03iLy one set of 
32 locations then the complete DPT has to be over in 5 ps. 

In this implementation it is proposed not to use any 
hardware multiplier. Instead a table look-up with a adder 
and hardwired shift is proposed to be used, Por each weight, 

16 in this case, we have a table with 16 entries. Each entry 
corresponds to the product of ten weight with the integers 
0 through 15. The integers 0 through 15 represent the 
numbers that can be represexited by 4 bits. If we have a 12 bit 
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numtier in whose product with the weight we are interested, 
instead of doing 11 shifts and adds (assuming that the weight 
does not have less than 12 bits) we take the first 4 bits and 
use it as address for the table look-up. The table contents 
are read out into one of the inputs of the adder. The eight 
most Significant bits of the output goes to a latch which is 
connected back to the other input. H'ow’ the next 4-bits are 
read and similarly acted upon. After last four bits are 
acted upon, the adder output goes to an output register (taking 
care of the sign bit) tne latch is cleared. Now we require oust 
3 adds. By making 3 copies of a table and dedicating each to 
a 4-bit portion the address look-up can be simultaneous. The 
output being held in registers/latches which are switched to the 
input of the adder in sequence. 

The flow-chart for the implementation of this DFT as part 
of its simulation is shown in Appendix A. 

A word about the DPT simulator. Its purpose is to deter- 
mine what effects different register size, memory (table) word 
size and quantization levels will have on the result as 
compared to the results from a DPT using a software PFT based 
on the same algorithm. Tovrards this end the parameters DATPOS 

rosMtn posncM 

and MEMDO B ( jIBiyiPes has to position tne memory word in the 
most significant bits of the register) are used so tiiat a 
smaller memory or dataword may be appropriately positioned 
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in the register word to get the best result, Fixed point 
arithmetic is used. It was found that by scaling the data 
word such that maximujii data word was less than maximuiii register 
content by a factor of 16. i.e. MTPOS positioned the data word 
in the register after leaving the most significant four bits, 
no overflows resulted with no reqLuirement of scaling in the 
middle of the processing. 

¥e now consider the timing for addressing the table and 
operation of the adders which will determine the memory require- 
ment and the number of adders. 

The schematic for the multiplier implementation is shown in 

J'lg. 

From the spec, sheets of I)M74S271 (256x8 bit ROM), SN 74/8283 
(4 bit adder) and SN 748412 (multimode- buffer latch) talcing 
maximum times, we see that 60 ns are required for the table 
look-up then 27 ns for the latch output then 22.5 ns for 12 bit 
add again 27 ns in the latch and then 22.5 ns for adding. Then 
the output from the output latch is taken after 27 ns. This total 
comes to 18S ns, further 13 bit addition we take 25 ns. This 
makes the total 214 n.s. If we take 30 ns for each complementing 
we get a total of 274 ns. 

Now of the 16 weights two, namely for 0° and 90®, can be 
hardwired, since for 0® the weight is 1 and for 90° the real 
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and imaginary parts are swxtclied (after changing the sign of 
the real part). Thus we have 14 weights which have to be 
multiplied and 5 (is to do it in, thab is for each weight we 
have 5 (is/14 550 ns. ¥e therefore have a leeway of 

(550 - 274 ns ' 55 ns). Thus it ifould seem possible to 

multiplex the same for all the 14 operation to complete the 
first stage of the DPT. The second stage will require 12 such 
multiplications. The third stage v^ould require 8, The fourth 
and fifth stages are hardwired and would require only as much 
time as required for doing the bubterflies. 

Thus for the first three stages we require three multipliers 
and since the third mulbiplier is required only for 8 opera.tions. 
It may share the load of the 1st stage mulbiplier. Thus we will 
have two multipliers doing 11 weighting operations and one multi- 
plier doing 12. Thus for each complex weighting we have 
approximabely 400 ns ( 5(is/l2). Thus we have a leeway of 

400 - 275 = 125 ns, for each weighting operation, to take care 
of the multiplexing and timing. 

A list of the o hi - p s required for the above implementation 
of the DTT is given at Appendix B. 
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4 THE MTI Filter 

A moving target will produce a doppler frequency which 
■because of the sampling at the p.r.f , rate will cause an aliasing 
about p.r.f. The blind speed for the system considered is 

* p", r ' . f , ~ where !^ = 2 m,^ p.r.f. = 400 pps 

^R ~ ^ 400 = 400 m/s = 1440 km/hr. 

It is seen that the first blind speed is very high caused 
by the use of a v.h.f . frequency. This is an advantage. 

To detect a moving target and suppress the response from 
stationary/low doppler velocity targets the consecutive outputs 
obtained at each of the 32 angular positions, is passed through 
a digital filter. 

The flow chan for the digital filter is shown in Fig, 

It IS based on a 3 sbage casca.de of 2nd order sections in 
direct form ^ as mentioned on page 45? Ref. 22]. To determine 
the values of the various coeff icienfcs the following method was 
adopted, A unit sequence (in time domain) i*c» 1 followed by 
255 zeros, was given as input to the digital filter routine. 

The output sequence from the digital filter routine was fed 
to a FFT subroutine. The magnitude of the 256 point frequency 
domain output represents the filter transfer function at each 
of the 256 points between 0 and p.r.f. The magnitude is fed 
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Filter 
Step 1 

Step 2 


Step 5 
Step 4 

Step 5 

Step 6 
Step 7 


Flow Chart 

,* Read in values of coefficients IRJia, irai42, NUM3, 
DNOMl, FH0M2, DR0M3. 

: Read in inpub sequence into stage 1. 

Ooiament ; The routine uses enough memory so that at 
each stage the complete sequence is retained instead 
of the previous two values only. It could just as 
easily have been done with two storage location. 

But for purposes of debugging and the initial 
conceptualization for coding the complete sequence 
was operated and shifted down. Thus IRTER (7,256) 
stores the values at bhe 7 stages shown, 

: Process the first pulse through the 7 stages. 

: Process the second pulse through the 7 stages. 

The delay units have now been filled, 

: Send the subsequent pulses through the 7 stages 
one by one. 

: Delete or weight :he first 3/4 oilputs. 

: Return, 


9Cb) 
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to a plot routine which displays the log oi magnitude versus 
the discrete frequency on the CRT terminal of the computer. 

The method is analogous to tuning a hardware circuit 
using a signal/pulse generator and an oscilloscope. Heedless 
to say it requires a little patience. The coefficients of the 
numerator decide the cut off region by appropriately weighting 
the denominator a pole can be positioned so that the transition 
region slope can be made steep. The frequency response of the 
filter and the coefficient used are given in Appendix G, 

If the MTI was required over a range of say 100 hms 
involving approximately 100 range cells each range cell 
would require 6 storage locations for the digital filter stage 
plus 1 for the input. The digital filter would have one PRT 
2500 [is) to attend to all the 100 range cells, i.e. 
approximately 25 ps per range cell. There can be 32 such 
filters to attend to the 32 angular positions. The MTI storage 
location required would be 32x7x100 = 22,400 storage locations 
for complex numbers. As bhe DFT unit puts out a 26 bit 
complex number at each angle position bhc total storage 
requirement would be 22,400x26 bits. 

By using a look-up multiplier as outlined for the DPT 
or a hardwired multiplier it should be possible to do each 
26 (complex) x 8/9 bit multiplication in 500 as.-~for the 
8 complex multiplications this means 6 ps. As each cell has 
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got 25 }is to 136 processed there is sufxicieiit time to complete 
the 12 complex additions, 

ihe MTI implementation would be the same as a digital 
IITI in a conventional radar only in this case 52 of them 
would be required to operate on the 32 angle outputs. 

The output from the last adder is in complex form^this is 
to be converted into a magnitude. There is 25 iis to do this 
using tvo 16x16 bit multipliers and one adder it should be 
possible to do the two multiplications (squaring) and additions 
in 1-2 [IS, ■ , . 

After checking agadnist a threshold this information 

IS fed to the 8 sectors xirhich are aliased. A logic circuit 
which determines the sector illuminated by the transmitter 
based upon its synchro information will enable the correspond- 
ing sector and disable the others. The details of extraction 
of angle information and the logic circuit are not discussed. 

An outline sketch of an HTI filcor has been given. The 
design has not been simulated, though the algoritiim has. It 
IS felt that ground returns at ranges in excess of 30 kms will 
not be excessive axid could be taken cal's of by using STC 
(sensitivity time control) and biasing using a clutter map. 
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GHA-PTBR 4 
SBIUMTIOB 


C 

The program for the simulation is given at Appendix &, 
The algorithm is given below. 


Step 1 


Step 2 
Step 3 

Step 4 


Input antenna parameters. Input DFT simulation 
parameters. Target and clutter parameters. Switches 
LIMTE2, IIMTR, TAPER, GOW. 

Load the R0I4s for passing lo DPT routine. 

Form the angular information bo be passed to the 
PLOT sub-routine, 

Perm the parameter to be passed to the RXROUT sub- 
routine which simulates the limiting action of the 
IP amplifier. 


LOOP : 
Step 5 

Step 6 
Step 7 
Step 8 
Step 9 


Gall noise subroutine (returns value for each of 
the 32 antennas). Add doppler phase shift. 

Veebor addition of clutter target signals and noise. 
If LHiTR call RXROUT. 

If 1IMTR2 call RXRl. 

antenna 

If TAPER then begi n if/unserviceable ¥T = 0 else 

WT = 0.54 + 0.46 cos end. TAPER the 32 

samples. (nNT gives the antenna position- it will take 
values -16 to -t-16). 
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Step 10 : 
Step 11 : 
Step 12 ; 

Step 13 : 

Step 14 : 
Step 15 ; 
Step 16 ; 
Step 17 ; 


Store xn array. 

If all antennas not done go to Stop S, 
then 

If CONY/begin call A2D roubine, call IPISH'I end. 
Else call ET routxnc. 

Determine the amplitudes and phase of the spectrum 
return by the DI’T routine. 

Store these Tulues, 

All return processed? If not, go to LOOP 
Plot amplitude. 

Plot phase. 


Step 18 
Step 19 


Step 20 
St ep 21 


Call DIOEIL. 

Determine magnitude. Average the magnitudes over the 
complete sequence. 

Plot magnitude aga.inst angle. 

STOP. 


It IS assumed that the extreme left antenna is the phase 
reference point and thax targets and clutter are all located 
at equi delay locus from this ref ^r once point. This would 
not in any way effect the results of the simulation. 

¥e can simulate the effects of noise level, local osci- 
llator instability, inaccuracj^ m antenna position as also 
the DPT parameters. 
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The Sinrulaiiion is only for one set of targets in a 
given range cell as the simulation is for verifying the angular 
resolution system proposed, 

.1 The Foise Subroutine ( Appendix DD) 

White Gaussian noise samples should have a constant 

magnitude in discrete freq.uency spectra. Therefore we take 

such a spectra. The phase of the spectral components we take 

to be uniformly distributed from 0 to 2%, This complex spectrum 

we now inverse transform to obtain the tune series. This is now 

our noise se<iuence. That the real and imagma 3 ry components 

will be Gaussian follows from the central limit theorem, the numbej 
of spectral components taken being 32 , 

4,2RXR1 Subroutine ; 

This subroutine simulates the limiting action of the IF 
amplifier and the coherent detection by multip3.ying it with 
cos (wt) and — sin (wfc). The 1000 points are taken over a 
cycle multiplied and summed to determine the output of the 
coherent detect on The symmetries in the sin and cos functions 
exploited as can be seen from the routine given in Appendix 

B. 

4.3 EXROUT Subroutine : (Appendix F) ; 

Utilizes a 4th order Eunga Kutta integration to do the 
same 30 b as RIRl. It ms used to check the performance of 
RXEl, The time taken by this is about secs whereas RXEl 
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takes about dd. secs. 

Signal to Noise Ratio ; 

¥e must determine what must be the spectral amplitude 
of the discrete noise samples that must be added to the signal. 
These discrete noise sample are taken at the PHF rate that being 
the rate a,t which each range cell is sampled. 

Considering a target of Im® at a range of 150 km. The 
transmitted peak power is 750 kw. If the transmitter antenna 
beamwidth in the horizontal plane is 15^ and the vertical 60°, 
Then, as 

Q 1 g— [Ref. Slcolnik Radar Handbook 


S “ 15x60 - ^5.5 h 16.6 db 


The signal peak power returned from the target 


_ 7 50x10 s X 45^5 x 1 _ 54000x10^x10 

If the receiver dipole has a gam of 0 db and an efficiency of 

0,5 its effective capture area = 4 G«p 

AtC i / = 2m 


X 1x0.5 = 0.159 m' 


Therefore the signal received at one antenna = 0.159x43.3x10"' ' wa" 

= e.SxlO"^"^ 


watts 
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¥e have the relation between output signal to noise ratio 
and the input signal to noise ratio as 




where is the noise figure of the receiver by 
definition. 


(S/IT}^ = 


mm 


0 n 


where B is the bandwidth of the receiver 
k ~ Boltzmann constant 

kT^ = 4x10“^^ watts/Hz assuming = 290^K 


Assuming = 9 db and since pulse width is 

T = 5 tis B I = 400 kcs, 

0 


¥e liave 

(S/¥)^ 


6 . 8x10'"^'^ 
4xl0''^z4xl0^x8 


= .053x10"^ = .005 

Therefore, E/s ~ 200 

If = 200 x S 


Since the manner in which we generate the noise leads to 
32 spectral lines in the range 0 - 400 c/s (prf), each of 
these spectral lines will have a power of 200xS/32 = 6.25 x S, 
in other words their amplitudes will be 'f§.25xainpl, of signal 
= 2.5 X Sig, amp. 



50 


Thus the spectral amplitude/eigiial amplitude ratio must he 
^ 2,5 when tney are added prior to simulating the limit IP 
amplifier and coherent detection. 

If only 16 pulses were processed then ratio 2.5xV2 'V 2,85 « 
If 8 pulses then it snould he 2,5x2 = 5. However higher 

ratios are considered to partially offset the fact that target 
seintillation has not been taken into consideration. 

,5 Clutter to Signal Eatio ; 

As the transmitted beam width is 15° ( l/4 radian) the 

patch of ground at range R illuEiinated during a. pulse width, 
assuming the complete beam hits the ground; is 

£■10/4x2 = Ss6x ig:k;ooxio ^ 

= 0,225xR sq.km. 

at 100 kms the patch would be 22,5 sq. km. Taking the 
scattering cross-section as -50 db [24]^ llie clubber cross- 

5 

Section would therefore be 22,5x10 sq.m. Thus compared to a 
Im® target at the same range the total clubber power would be 
22.5x10^ times greater or the total clutter voltage v;ould be 
150 times greater. The simulation has provision for 5 point- 
clutter whose aaimuth separation from the ta,rget, relative 
strength and doppler can be read in. Thus each clutter point 
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will have to have a voltage approximcely J x 10^ ^ 67«5 tii 
times that of the signal voltage. Admittedly this is much too 
simplistic a clutter model. Hovrever, it is unlikely that at 
the ranges under consideration the complete beamwidth would 
illuminate a continuous patch of clutter. The above however, 
gives an idea of the clutter signals to be used, 

4*6 dumber of returns from a Target 

Since the transmitter antenna rotates at 6 rpm and has a 
beamwidth of 15° the number of returns from a point target for 
400 transmitted pulse/sec, 

= = i|22 = 166 pulse. 

The simulation processes 8/l6 returns this is to facilitate 
the assumption that over a small section of the beam the beam- 
shape and scanning effects caii be ncolocted. 
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CHAPTER 5 

CONCLUSION AND SUGGESTIONS 

5.1 CONCLUSIONS 

The simulation has brought out 

a) The requirement of preventing the signal amplification 
stages from limiting. 

b) That under expected noise levels and clutter levels (30 dl) 
it is possible to determine the azimuth angle, 

c) That the full conversion range of the A-D converter should 
be utilized by incorporating appropriate biasing, 

d) That register word lengths of 12 bits (11+1 sign), and 

8 bits in memory tables should give adequate performance, 

e) Owing to endfire resolution being poor two linear arrays 
each of 32 elements and mutually perpendicular would be 
required. 

In conclusion let us list the major system features as 
envisaged , 

a) Two mutually perpendicular antenna systems so as to 
avoid loss of resolution along the endfire direction, 

b) Synchro information from the tiansmitting site to be 
transmitted to the receiving site for deciding the 
display sector. The other information is trigger and a 
lock pulse for the coherent oscillator. 
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c) RP amplifier to be antema mounted. 

d) An lAG-C circuit and/or a biasing voltage proportional 

to the clutter bo prevent the IP amplifier from limiting. 

e) Coherent detection followed by BIT unit and possibly 

a HTI unit followed by display unit controlled by logic 
circuit for determining the sector, Ihe max delay caused 
by the DFT and MTI processing is expected to be 50 jis, 

f) A switching arrangement which will enable the mixer and 
subsequent stages to be shared between the two mutually 
perpendicular arrays. 

An overall system block schematic is shown in Pig./o 

An apparent disadvantage may be that for given transmitted 
power the ret^mn from a target will be less for the antenim 
with a wider beam (lower gain) while the groiund clutter 
assuming it bo be uniform would return the same power. Thus 
the input signal/ clutter ratio would be poor and hence sub- 
clutter visibility at the output of an MTI system would also be 
poor. But what has been proposed is for a system where the 
beam has been made as narrow as operationally feasible and one 
IS now trying to look inside the beam. The input signal/ club ter 
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ratio tiius -will in aiiy case te tiio boat possible under the 
circumstance. 

The scheme though proposed for a vhf system could in 
principle be applied to system at higher froluency, lor example 
for a radar m the L band ( ^ = 0.3 m) wiLh a beam width of 2° 

woixld require sampling at approximately .07^^ samples/motre 
If sampling frequency wore divided into 8 disci'ete points the 
8 spatial samples would be epa spread ouo over 8 x = 240/7 ' 34 
metres. At 60° from broad side such a sysrem would still give a 
beam width of approximeitely 1°. At broad sido the beam width 
would be approximately 0.5°. Of course xhe number of sectors 
which have to be resolved has increased, 

5 . 2 Suggestions 

Since in any array there will be some interelement 
coupling one convenient way of compensating for this would be 
by multiplying each of the samples by a compensating factor 
before it is fed to the DlT.unit. 

By taking the concept of spatial frequency and the 
band limiting action of the transmitter antenna, we have 
by applying conventional signal processing, namely a been 
able to obtain angle informtion from within the beam. It 
may be interesting to consider other signal processing methods 
(e.g. an elliptic filter) which will give non'’rox7sr passbands 
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and could be used in a given direction. 

By introducing a threshold check beioro A“D conversion 
it may bo possible to get more tinie for processing by ro] ecting 
those sample sets vnere a given nouber of samples fall below the 
threshold. 

It may be worthwhile to implement bhe DPT as well as the Mil 
filter using CCD devices. This would at least cut down on the 
A/D conversion. Electro-optical array processors have also been 
reported [ 25 ]. 

An important aspect thafc will have bo be considered in any 
final design is the monitoring of various parameters and the 
incorporation of automatic self- check features. 

Evaluation of probability of detecuion and probability 
of false alarm has not been done. The simulator that has been 
coded may be extended to determine these for different values 
of SiTR, under conditions of sciiitillati-ig' target and a more 
approprn te clutter model. 

5,3 Discussion 

How does this method compare with tho one proposed by 
DEN DAVIES et al [25]*^ Ihe method proposed by Davies et al 
uses a swept frequency local oscillator to introduce the 
phase, the mixer outputs are then compared in an IP beam forming 
network to obtain the received beam. As the transmitted beam 
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rotates the local oscillator phase shift uould have to be 
adjusted so that tne received beam seuiiigs within the trans- 
mitted beam. Secondly assuming that 32 received beams are 
formed in both cases. In the case where beam forming is done 
by sweeping the local oscillator only l/32 of the total energy 
(over one pulse width) received from a target is utilized - 
assuming that the received beam looks in each of the 32 position 
for the same length of time. Thus the video output after detection 
and lowpass filtering will be much lower. Whereas in the proposed 
system as the total received energy is stored and then appropria- 
tely combined using the DfT no such loss as caused by ' distructive’ 
combination can be expected. Thus to get the same performance 
of detection the power bransmitted in the case of the former 
method would have to be 32 times greater or the max. range 
reduced to 1/ *^3"2 0,42, Using the method proposed of course 

would involve coherent dotoccion for oaC'i of the 32 antenna out- 
puts which would not be the case with the method proposed by 
Davies et al. 

[27] proposes a holographic radar suitable for high range 
resolution at small ranges, the application being measurement of 
ice tnickness. The radar operates in 0\I mode and uses the 
holographic effect for determining the range (depth). The present 
proposed processing resembles it in that the output after coherent 
detection is converted and stored and subsequently the beam is 
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FLOW CHART OF DPT SIMULATOR 

Step 1 Input the following parameters : 1) Ho. of aerials/ 

spatial samples (lORD), 2) HSTAG- log of spatial samples 
in base 2, J>) Ho, of quantization levels (KWAHT) , 

4) WRLSIZ wordsize corresponds to the number of bits in 
the registers whiC5h will hold the output of each stage 
as also the adders, 5) Ho, of bits for each entry in the 
table MEMl^EL, 6) POSIfflM appropriately positions the bits 
of the table look up in the registers, 7) DATPOS positions 
the data bits m the register, 8) DCDBIT the number of 
bits used for address decoding. 

Step 2 Load the tables with 2 3E s DCDBIT entries for each of 

IORD/2 complex weights for purpose of simulation rotate 
the bits corresponding to ICBMPOS, 

Hote ; When calling from main routine Step 1 and 2 are 
passed as parameter Step 5 is not used. 

Step 3 Set random number generator generate (lORD) real and 

imaginary samples and quantize them put them in arrays 
STOR and STOI. Retain the unquantized numbers in array 
XTEST for verification. For purposes of simulation the 
contents of STOR and STOI are rotated to correspond to 


DATPOS. 
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Step 4 Commenoe outer loop of PET sxmulator. ISTAG- taking 
values 1 to NSTA&. Kl = ISIAG- -1; K2 = 2 3 e * II,- 
I = 2 ae 3E {ISIAG -l)/ir2. H2 is now the number of loop$ 
in each stage of the ISIaGs stages. I is the number of 
butterflies in each of the 12 sections. 

Step 5 Do the loops corresponding to each stage. lo, of loops 
for n taking values II to 12. N5 = ( H - 1) * N * 2. 

13 will give the address for register for the first input 
to the first butterfly in each of the 12 loops. 

Step 6 The inner most loop. For 12 taking values 1 to N - 
I is the number of butterflies as given above. 

14 = 13 + 12, 15 = 14 + I, II = 2 * * (ISTAG - 1) i 

(12 - 1) + 1. IN IS used in determining the appropriate 
weight in Step 

Step 7 Real and imaginary butterfly operation. Real 
BTP1YR(I2,1) = ST0R(I4) + ST0R(I5) 

BTP1YR(I2,2) = ST0R(I4) ~ ST0R(I5) 

Imaginary 

Ba]FIYI(I2,l) = ST0I(N4) + ST0I(I5) 

BTPLY(I2,2) = ST0I(I4) - ST0I(N5) 

Step 8 Obtain the sign of BTPLYR(I2,2) and BTPIYI(I2,2) as 

these will be involved in multiplication by the wei^ts. 
Step 9 Check for error. If any of the butterfly outputs will 
result in register overflow then go to ERROR, 
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step 10 Put the sum output of the butterflies back in 
S!rOR(N4) and S!r0I{N4) respectively. 

Step 11 If IS^IAG eq.ual to NSTAG then store the difference 
outputs also as no weighting required in last stage 
and go to Step 24, 

Comment : ¥e have 4 adders two each for real and imaginary 
parts. They each have two inputs represented as ADIRll, 
ADIR21 for Ho. 1 real adder. ADIE12, ADIR22 for Ho. 2 
real adder, aDIIII, AD I I 21 for Ho.l imaginary adder, 
ADII12, ADII22 for Ho, 2 imaginary adder. The corresponding 
outputs are ADDRl, AD0R2 and ADOIl, AD0I2, No,l real 
and imaginary adders handle the product of the real part 
with the complex weight. Ho. 2 real and imaginary adders 
handle the product of the imaginary part with the complex 
weight , 

Step 12 Clear the inputs to the adders, 

LOOP 

Step 13 Cet the output of the adders. 

Step 14 Shift each output to the right by the number of decode 
bits DCDBIT, 

Step 15 Extract the least significant decode bits (which have 

as yet not been used) from the difference output of the 
butterflies, i.o, BTPLYR (12, 2) and BTrLYI( 12, 2) . 

Step 16 Use these bits as addresses for looking up the table 
corresponding to the appropriate weight. 
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Step 17 


Step 18 

Step 19 

Step 20 

Step 21 
Step 22 


The res-olt of the table look up is put into input 1 
of the appropriate adder. aDIPJLI gets the product 
of the real part with the real part of the complex 
weight, ADIIll gets the product of the real part 
with the imaginary part of the complex wiweight. 
Similarly, ADIB12 gets the product of the imaginary 
part with the real part of the complex weight, 

ADII12 gets the product of the imaginary part with the 
imaginary part of the complex weight. 

If the most significant decode bits have not been used 
go to step 13. 

Form the adder outputs and affix the appropriate signs 
saved in Step 8, 

Form the real and imaginary part of the product by 
appropriately adding and subtracting the four products 
obtained as output of the adders in Step 19. 

Check for overflows. If ifany overflow go to ERROR. 
Put the result in appropriate storage, i.o, ST0R(E'5) 
and ST0I(N5). 


Step 23 END of inner loop. End of middle loop. End of 
outer loop. 

Step 24 Do bit reversal and shuffle the storage accordingly. 
Step 25 Scale output down by shifting by DATPOS bits. 



step 26 
Step 27 
Step 28 
Stop 29 


Form the magnitudo of the complex output. 

Call FFT routine using input as ZTEST for verification. 
Chock error of corresponding samples. 

Stop. 
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MflDWARE REUUIEED 
POP lEJLTIPLIER 

It em Qty 

ROM DPr74-S271 12 

8 bit latches SN74S412 6 

12 hit latches (two SN74S412) 12 

12 hit adder (three SN741S283 4 hit adders) 6 

15 hit 2' s conplement er 12 

1 hit latch 2 

SIGN logic circuit 

12x1 multiplexer 26 

Output demultiplexer 26 

Butterflies 

Adders 5 

Subtract or 5 

16x1 Input multiplexer 5x2x2 20 

2x16 Output multiplexer 5x2x2_ 20 

Stoiage/latch 

12 bit latches 2x52x6 584 

Clock and timing circuits 
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1 ( 3 , = ' ,13, ';>A n-'iis-' , Ti ) 

hi nip I i ! i i n i I i n u i 


^ 1 3 


IORD=' , 


L 'IH’K ♦, (riiiTuCK) ,('.rn C<V,) , tiJHO) 

CAl.b FTtXi'EGlMSlAG 
CC 1ii'K * , (XVEST ( n , 1 = 1 , LfJuOJ 
C fJG aft I=1,I.URt; 

Cbb ()=CAH5CXT-jSr(I)3 

CC rVPK ’►,CXi.AuCD,fMGCU,I,I = l,inKD) 

C CAGi, FTESi ( KT'EST, gGTAG, ,iaHD3 
c SvPPVKso, ?Do sa 1=1, U 

C58 , KLi:J'''.r(I J=0 

C or; b'l r = i. ,IURU 
c AT .\i;i.r)=rA4SCX'XKS'i({,)} 

C .■:nR=(X‘'AC:( r. )-FLUA'r(MAUC X) ) )/XMAG( I) 

C n%P M:=iK IX CABS (ERR )* 100. ) 

C if ( IPKH.CG.iO JGO m 54 
C f^uOtTtlERH+ilsKilUOTdERR + D+l 
0 Gu 10 57 

C 59 KDII.'JT (11) =K0IINX ( 1 1 ) + 1 

C SCREWHsSgHE^^H + (ERk)4<'2 
C57 Cat'H’lWlIE . 

e SlGERH«SORT(SOBEFa/KL.jATCiOHDJ) 

C type 802, SiG£PR,TUPO,(t, KOURTCI), 1=1, lU 

C802 FOR/ai C.lX,'HriS UEV I ATI OMs ' , FI 0 , 5 , 'NO , OF POINTS IN THE FFT',13, 

C U/lX,'Nn, OF PJI.JIG iM Si 3, 'IS', 1 3)) 

C CYPK T,(XMAGC l),MAG( U,I,I=1,IQRD) 

Opl.l, SA rfRA'HISA^) 

C 1 «f-f. as , USfiV 
i IS AV 

ro>jiiGUE ■ 

■O") r<;rMAI( IX, 'XSAV ='I12J 

C.ju I l.a/f;; 

/' ;■ i !•' . ■ 

1 < , 1.1 ' ■ ■ 



-B) tr- 

' ' ' 1 .'-.t KTKSrfX, IN,M) 

IIKdO) 

C<).'’k! K.X XtN ) ,w ,lvrFLlfl ,BrFLYii,XEH,TST 

'i ' . ■ P 1 s V , ♦ i , 1 '» 1 !> 9 ? W’ = C ^ P t , X ( C ,0 5 C T U P I / F L Q A T ( r'i ) ) , - S 1 M C T U P I / F L Q A T C N) ) ) 
‘"y=:H/^.V'JSrAG5=IiV 
i'-!,. / i7=t ,«.‘3-,rAG 
I M iP-1 )JI)n 6 I irl , lY 
u,='-x/i ) ♦uti ? izzsizK I x-n 

- i).j ’•i ii=3zazz. 

1 11 = ] i 

1 i=n t + rx 

C type »,'<( 11 ) ) , XlU J , 111 ,Ii 

B iKi.y isx( n i )■( xc u 1 ; !u’FLY2=X( 111 )-x (ij ) 

1\S1=V. *♦( (I) i-J Z J 

/-) A t U 1 )=;-iTM.V 1 ; XU l) = BiFhy2*TS'r 

1 ) 1 B 

1 Y-) Cur TTNUF 

‘,1 T'.-'T! ■:'!>£ 

n f , ' i f p;!'/ 

I ( . iUPji-' 

'■ T3 Tl-u. 
i 1,4 = 5 , 

n V I U) ) = 1 ' 

'■ j = l, 7 1B = T1-1 

II 1 rRM = i 5/27 1F( rTt,*'','vy,0) GU TO 10 

1 ix( ^'n=(’it'u I'h ,'/>;= i~“i +1 '* ’9=iTK^?r.ij to u 
10 nxu n=ib 

K = 1 ■ , 

i ! 12 F I = 1 , N .S r A ; M 1 = K I - 1 7 *P,f P M = 0 S 1 A G - K 1 1 7 F = K + 1 1 X c N T E M ) * 2 * 1= C K 1 1) 

12 Cut I I ',10-; 

C '* » !''■ S'. 'U » I 1 

. ( . u r i ? p I i .1 
' c-/ u > ; a( .u=:u n 1; X ( 1 n=TE'i ' 

, :• I V 

■ I ' ' » i 7 

‘ • u ‘ • U I , t = I , 



DD !■ 


f 



oi ..P.IIUTIHK i^OlSliifXY, 

Ur Ah Kii? ,X) 
j.iTLGKH M/iCn 
rnMP!,ia XV (W) 
c’ (’u-tpl.f’x 

i' ''i=X 2 ;lN -5 . 

C OU AO 1^51^32 

■ UO 5 I = l,i'i?FIsRANCi)UMj*i'iJpI; 

.XHI)=CMPjbXtCijSa*'l ) ,SlNtFI) J 
5 CONTTNUt: 

Call FT(xuiN,N) 

r)U 10 1 = 1 ,H; XYM )=Ci.»'^0GCrX(l))/FLOArCN) 

1 ■> Cbr' TTuot’: 

i’LRiv 

f.c H1Z.)=XYCi?>) 

(' ‘n: .'n'/l = l ,'(,*xYtTZt ) = TFM(lZn 

/ '! i 1- i' i'i 

t'c ‘--I,: , iXY ( i. ) ,1 = 1 , u 

C CAi b 1-"UX / , LV,R), . . 

c nv';=o, . , 

c i:-i (5 1 =1 , a;x (I , n=CAiiS(XY(i'))j y (j , 1) = I 
t AVC=AVGtXC I a ) 

15 CONTIWilt 

C AvGsAyG/FLUATCJ^) 

C Sl&«A=0.?l)O 25 17=1 /NaSlGMAsSIGMA+f/Cn.n-AVGJ^I^^Z 
25 CUNTJNUF 

c L lG,'lA=SORi'(.SlGMA)/F LnAT(!'J) M’YPF. ♦,A,VG, SIGMA 
C ' I 1 =!'i ; S' iMT= . FaIiSF . » M= 1 
' t /' Lb F IjjT k, S'JF' i f > f if A f A; A f A 1 

(■ -i : :r^ V 



I 



oHiMjT«’UK tp:sti ’ : . 

/SL LftiBOA ' ' 

'.'ifUCAl- MrrJ'rti? 

I11S9 

Tit'Pn; h 

;ATC jl X, 'It'ipn t lambda , phase, CUTOFF, Z.M' ) 

ACCKPI XI N, LAMBDA, PH ASE,CUTOFF,ZM 

CALL HXKOtiTlXJM ,f<AMHDA, PHASE, UirrR,nUTl ,CnTOFF) 

1 IMTP'^s.THrrE, 

)1.»'K102 4? ;)W1 =f4W/4;i)f.;L=TUPt/FLOAT{N.O ; N«2 = Mi^/2j NW3 = NW2 + L5nU 
FPKQr JOO. *10 ,>«'*'3/bAf!B!)A; n 1= L, / (FREO^FLOATCNW) ) , 

P'-'Ff'AsTBr I ♦KPEO 
1 t = ! -t 1 

' it; i=i ;;U=:,,>lu(FLUAi'C l-l) + PET,) jCOrCUSCFLOATCI^l )*DEL) 

!. 'I i=i:'-!‘i„.\a'!i,si ) 
f '* ?>: 

lift r»T .‘7 KPMJ, «ym fWw,MWl,l)FL,»yW2,NW3,MWlt ,xin, omega, phase, DT,W 

) , ri' , ,'/,M 1 

"\'-i 4 ,..t ■ V j .kK,A IM 

r' jH 

‘I’V kKHI. (-.a , gw I ,DFlj,NW2,hW3,MV/ll ,XIN, omega, PHASE,nT,W,RE, 
1 Ar':,CUT('-FF, Jl'’’TR) 

HKAL X(1024J ,XX(tL241 
HEAL J1TTH . 

COMP!,, EX ¥(1024) 

COMPT.EX WCNW1 j,i’EMR,Xj((li)24) 
tEX=tO 

t4J SO ■ 

r> Vf?.'xpi(A;jFf J,rTrP + TiJP( . 

' i t-' ' , ' 

' I".' , ■■ 'j r I' --'vsFi.uA r c 

^ ^ 4 At ~t i ’ ' -I i* 

I > A I 1 ( , , ' . ; I 

' A- ^^, >, i: i 1 

i t . f 1 ) J . vr.t’i FF! a( n=rilT()FF>t!3lGN(1 . ,X (I) ) 

•.’.,,.(11 

I i . .'A . f A r ; > A-' . ) ■ ■ 

< ' I 

*' I : *1 , I ’ < I ^ - I , . ) 

^ A 7- 1 • if.-# ,1 ? ^ C y « # # .1 

* A , » ■■ 

' ’ . I *' 1 i I I 

- i - t M I I - * r ! f IU( vi(.l ) } + 



V 



.(E 2* 

i'f'i 1 !i 

.i . 'K.'.CiiXU )“X(,4W2 + 1 )) + C0. ,-l .)»(X(NWlfXl-XtNW3 + I))+Tf:HH 

* rUNTlfinfc; 

«P,KAl,(,T'}'r«^H); A IMsAIMAGCTEMRl' 

P'-,Ki\’FVK|,nAT(fv«il) ♦?. ? A tMaftlM/FLOATC NW) *2. 

^ !)*» (’AT, I, KTa,1U,1024) 

C TYPt:; ♦■,¥(2) 

Hb'TliP'^ . ■ 

KKD 



i 



-F i" ■ 

.■ f l • I r 'IKST^ 

' 1 ' j 1 I. n Hi / A 

•■. t C?bf> ) ' 

1 . ^ 3. 14159 • 

1 s 

b Fuw'r f IX, 'I^PU'I, XIM , LAMBDA, PHASE, CUTOFF* ) 
ftCChPl lambda, PHASE, CUTOFF 

C f’Hi.J, «XWf'iJti>{I.nLAHbDA,PHASF,OUTR,OUTt,CUTOFF 3 
I. lf'VR? = .Tb.UE, 

/I «< = l <' 0 O,’i- O =N«i/ 9 ?DRD 5 TUPI/FL 0 ATtNW 3 ;Nw 2 =Ns«/ 2 ?NW 3 s!NF 2 + iMW 1 
pi: F. jiHU.*! lambda fDTsl ./(FRKO^FLOATCN^) j 

0 ,-!£(,A=Tur£»FHKC» 

I’. U s fDi 1 + 1 

LMi 10 ls 1 ,,-.Ml;RlsSlNCFLUATCD-l 3 »OFLJ ?CO=COSCFDOAT(l-t )*DEd) 

wu)»odfi xfc;o,s‘n 

10 COWTTNDF 

TF (ijf HTH? )CADb hXkH«w , o’Wl , DEL , 2 , S.W 3 , NWl 1 ,XrN ,OfPiEGA , PH ASF , DT , Im , HE 

1 , kK, AT? ,cni\:)FF j 

\ ’ i L I'"' i,.., p V.) 1, f t § ij ij f i ■ 

i I ■ . ■ , ^ 

SDi.l'ui* , l:..K RAFDU'If Xtb, LAMBDA, PHASE, OUTR,OUTI , CUTOFF) 
f.I'.Ml, 1 A irtOA, Xf 30 ,-l) , £( 4 ) ,LIMXT,DT,C 0 EFFC 223 

C I' ! * P i ' A 1 I 1 w ) , i ) , 1 I * 1 ^ U PI , 1 1 i C u C 

I t.n;' I- /A u ) 

Pi' Al, Xt) n ) , 1 D ( 1 j 

DODTCAL s* 

UAiA X/BOPu./ , 1 / 4^0 ./, i 9 / 4 !P 0 ,/,rOEFF/ 19 * 0 , , 1 2 * 0 ,/ 

c M=ino,’N=nHACi) = iou 

M »1 1 Ns) 

STAFTsO. 

Sw», FALSE, 

■ Kf-iK':' = 30 Q,»'if), ** 6 /LA’pBDA 
. = .,/*'’<Fv 

, f34KHFvi 

1 . i' 1= 1' 

; • . / 1 :• K’. 



I 1 i* .H' ) 

1 . , I , X.., m, :D'. , H) 

J-' I 1 , 1 ) , /T 



•f 2" 

j 

i I 1 5^2. /’i’ , i 

. . . i‘ F’ i f K « . 

W'lJ, PKfA, If ,lfN, LIMIT, COEFFjDT, START, XD»YD,H,(^,KNT1 ,S's<,OvE{;a, 
i t'-r ALK,XfiM, CUTOFF I 
. PM,i. f'LuTC .FALSE.,I ,XD,YO,MA,M) 

P liuTisYrCJ OO, 1 )*2./T 

P'lriEVljll )*2,/T . 

F-K'i’MHL ■ . 

M.ibFOti rii^K PK (X, Y, Y^,iaMrT,CiJEFF,DT,5TAHT,XD,¥0,MZ,NQ,KNtl ,Sw, 
ni2n-GH, PfiA.SFiXlN, CUTOFF) 

C H'UP'IH UHDCH RUMGF-KUTTA FUR DIFF. EQ, UPTO ORDER 20 

C Cl iflFF (Z? ) ARE (CLEFF ( I ) .D^^ao + COSFF ( 2) ,r*^i9 + . . » • • . . . + COEFF (20 j ,D+C0EFFC2n ) 
C +COF.fFl22) =FUNCCy,T) 

c o=uy/oi 

C X(l,l) ARE XHITIAL YALUES UF n=t»19 TO D 

c 3f (1) Initial value of i 

C XCi,n,Xa.n),X(I ,J) ,Yt2J,H3),Yt4j ARE INTERMEDIATE VALUES WHILE CALCULATin 
(’ 4'!,*' OPDF'K PU,-.;f;F~ia) ITA 

J H) 'vKf. vfli.UES r'F ,1**20 CALCULATE!^ 

^ ^ CIFJ 05 TriF*: INI IT. A L TIME 

v: s. i.-ir I;; Fii'>.HL I’lf'E. 

C I T I n V I T -CHEOF .',T 

U Af’L"N, ,,H ir.iT-) RE SFl’ lU TRUE IF APLUT IS WANTED ELSE TO FALSE 
C r. T-. •' if ijF Pul Fl'S HI.OTTEu AND IS THE DIMENSION OF XD,YD, 

C ..'Tf. . sf i,,T I n - oTaH n / (OT^KiiTll 
C AlT-Av'S THK RESULTS 

c (• Fui-crxn.', ourprocram rust »e provided 

WEAL CJEteTYx) ,rutT(22J 
HEaL X(2l , A) , Yt4) , 

C DlMEL’STuN XD(MZ,Ny) ,Yu(OZ,r-u) 

OIMENSIUH XUCl ) , SftKl ) 

n.IMEfJblOM SLOPE (4) 

1 '*( ' 2 0 ’ ' 0 1 1 ” t , 4 
( i 1 ) = U, 

Y I ! : 1 I = (' ( 1 1 J 
. < ' 'f I. ' O' 

■ .,i' .1 ! 1= i ,22;v-'nEi'(I1 j=COEFFMt ) . | 

■ t T’- ’ ' i' I Ri'. I 

, ; i '• 1 ) .1 ~ i » /. R 
! ! s: . 

! *' .if', Ti , r U 1 T'i 20 ■ 



I 1 ' t, T * 

, 

' ‘'i « 

, 1 ) Oil '*'0 30 

b 

™ 1 

XRMM 

0 ) 

3 

CUiC’ ( 1. i ) 


C’’ M 1, 


■; c 

. 1 y '1 GO ,10) 2n 


/, 1 1. , i 

) * X t; 1 r , M 

/ ' I 

<" = fM 1 





II 




! F C I 

\ , 

.bf- 

C, 22) on 2H 

i U 

!'l = lj4 2 





Ofi 


i<»2 rO 

ib 

OUET i.r i 

= 

COcO‘'<KM/CnF>-M ) 


{lOhf, i 


i , 

. /CUitFCi ) 


CUi-'F 

C i, 

) = l 

MO 


Ti(Pc: 

000, 

, (CaFFM)M=t MM 


TYPE il 

n FOHf^ATClH ,'i4,nMAb C0'*niT»04 FIRSt IK VlfO) NFX'f FROM RIGHT TO LEFT 
IT uf 0) ,n(UJ ♦4? — , f /I H 

5Uy FOMHATCtH .'CObl-' ARC t'J UKCRK A.5 I mG POWERS UF D RIGHT TO LEFT, HIGHTMOJ 

IHTMOS''’ rs crj.sr. •iR.xr ,5riu,2J) 

K*<=:0 

STr:P=KTAR’^ ■ , ' , 

. =0"’ I ' 

f,- f =<"• 

,1 y I I !•! 1' != ' ! , ' 

= STOP 

t ! I ' ' ' « ■ ' • t'' * ^ * 'j' ' TO i 2 J 

i\ si* iv ■!• ,l 

Xut! K,l, J3/l,iK . 

YDf f !■’ , 1 n J 
XLfFK)='nME;YU(R.K}=YU 1) 

123 CORTIPUE 

TFCSTEP .GT. LIi’AH') RKTUR.v 
1F(P! ,E0. 0) GU 'i’ll Aj. 

OU 40 1=1 ,0 ■ ■ 

•l'. =: C'TKFf r+ ri»X (.1 , 1 ) + TKH 1 

n , , ( 1 ) i-C.'Pv.'’ I'^FOwCt O'Ui. 'iA,RHASEf SW ,CUTUrF)-TFHrt'*C0EE(vH1. ) 

M- ( I f U2 } 

, I O { , ' / CM IK H. .4) 


bi.i".- ,1 'o; 1 1 ) 

!-'i M*1 

*5t. 

•c'^' i. - ^ ^ 
f 1 , 1 M' ^ 

{ C f. . ,i ) 

‘ . ! ' -X 1, !' # 1 1 

1 ^ - I f i I j ♦ . » i ’ M H 

.‘/,M 


F 4 


M i /I.* 

t . s; < , 

I- .fcO, U) GO TO 6 1 

wt’ Fu .t=.T,f4 ■ 

.f n/‘'.=Gi]5:K(X + l I^XCI ,2)FXFRM 

‘.1 V:U?J»C0wST*FU^C(TIMF,XIM, 0MFGA,PHASE.3W ,CUTUFF) -TKRrt-CneF Cxf+1 ) 

t 1 ( / )-ruFt (j+i ) 

rtW.tATil.i IS 'r4E14.4) 

5l,CPB')fiU2 ) 

!F(.'i ,ECi. 0) GO fU 71 
!'0 VU iBj ,H 

iCI,i1= SI.OP+Di/^, + K(r,1) 

7u ai.u»"=X (1,2) 

71 yU3)=Yl(U 4 8G0F*l)T/'>. 

IT'wwsd, 

TF('^ 0) GO TO 81 

UO Hu Isl 

no I’EKM a TFOGvi-H CuKf C I + 1 ) ♦X ( I , 3 } 

HI yw( :naCO.'’.S U, ,C(TiME,XIif, omega, phase, sw, CaXUFf3-TERM-COEF(J + l i 

1+Yl f i1-f’.jKF(.j + 2 ) 

,'A(I,4) IS ',4Et4,4) 

i' I , f' i; //.,+- T i ! 'i E 

SOI i' = v M’ 3 ) 

U-'C - .f.-:.', ■)) GO TH 4 1 

'U: 

,A ( I , . ■'a ■; . o L( , rj 

,5t!,0-s-''( f , 3 I \ 

'll ( 1 ( 4 Ja M 11 .1 4.sr,u>'»0r 

jp.t'''’ a'.», 

U'C?'* .GO. uJ Gu ro lut 
on 100 j=i ,/4 

100 TERM arFKM+COKFdl-U-XXt 1,4) 

101 y>U 4 )»COi<iST ♦FUi»,C ( TIME, A Ii< , OMEGA , PHASF. , 3W ,C UTUFF ) -TKRM-C0EF(G + 1 J 

1*Y1 f4)-ruEF(a + -(!) 

:i--H I -lOr. AT( 1H ,'A(I,4) IS *,4Ei4.4j 
I 1>.S !,al,4 

1 4; < f- 1!. ) a 14 (G,J 

i0.''C ,.)) no TO lit 

, ■ " li ' ’ = I , I 

■ K M'af o, . 'Gfi )+,Su'O'i«:r2)+2. + 2.i3LaPEC.l)4SG0PE(4))*l./6,l<Dl 

i » M i ) 

; . n I I I I,, S 1 , 4 

! ' •' ,0 ^ I - 1 = X ( 1 , i o 

r, ■ a I , n r. , 

n ,ii*i= f Cl J+PdOOUPFC^I+^.^SEOPECJJ+SUOFF.CAD/A.^DT +Y1C1) 



: 


»iF5- 

^ f’+OT ■ 

X*'n"s!K‘JT4t 
G'' !"■) 1?0 

MO 

['K?^T, FUNC'TTh'J' Fnf''rc TIMS' .XTK' . SW ,rnTnpi7) 

A, NflLKanf^FnAliTlf^F; ’ 

)( =X T ‘■'i 1"'0 S f A '^raT> F4 -P ha S F I 

rFCAPRCXI ,<|F. rilTnFFIXsrilTTF'S'^RTnM f 1 . .Xl 
IFCSO nn Tn 5 
S'l!MraX*CnSC A’^ni.EI 
RETURN 

<S riHXCaX^STNC hHCl F) 

RETiJRAI 
F M n 






G* I 


C .S t'i ■ t i'RSTil 

•m:; , ^ i),FHFy(256,l),SPMAG(256,l) 

■ ■ !■ >. vAii) 

HI..;-, X C <3 ) ■ 

' f':jRT 

~ 'i . 

i.J*- ; = J*1 4159*2. 

RK'i-l./ 4 O 0 . 

PmJ = '^U,*TuPI/360,0 

1 W K 8 ■ 

UU ip J=J,M 

RKF.;( j, I IsFLQATC J-1 J*400./FLDAT(>1} 

C Fhr-Fl.OATlJ-l )*6. 

I I' fui' Tl i'dJE 

Ih.SEOd , 1J = 1- „ • 

CU 51 *7, N 

C '! JM=KLnA'r(I,-l)*PPT 
e 1 rt.SE0tl f 1 C l'UFr+FHE*Ti(^' + PHJ) 
l-.5E0(I , 1 ) = 0- 
5 Can'lir.niF 

C T IPE 4 , U-vS-:OlK , 1 ) , K = l , n) 

CAT!, PTCrn L CT -kSEC ,S‘U 

(■■ . I, T .SK'..HK, n ,K = 1 . (O 

c 'i 'i'"-"-; * , 'li iM f 
C' l > 

1 ' .1 1 I i I 8E 

i. 1=1,' 

5r'Kv' ( n =r .f ux I j 'iFfcAH I. , n , 0 . 1 

lb C-'-. IT' 8 is 

c’^ii.u r- nr. PL c, IF,'") 

lu; 2(‘ isj,^ 

ftfMAfUl , J ,>sl0.*ALuG10(CABS(SPEC(l J ) ) 

20 CUNTToiUt 
MA{1)*756 
(’•lAf^T 9 Of! 

KfjRIVJ at C //? X "•"•“I. N DBS ’•-mmmmmmmimmmmmy ' ) 

('Ai,,, i' !Arr(SbHnM,FREO.SPMA,G,MA,ft) 


fj'.,. bl',,Kjb(!X-.Sr,0, n 

* tl n , JA A I n , .A-nT t 3 ) ,0 n 0.- n j) ,DnHn7 ( 3 J ,UMUM3 (3) ,K 

^ 1 . Ar.r; i, , 1 i 

> ' , ■"H, A (7, '/'•jH ,J 



UA ? 


'Ki- .i! . /,7i5J 

\ / ■ M L : ' ! / 3 ^ i I* / / * ' ' ' / o (i ^ 1 , / ^ U * / ^ j j 1^'’^ 3 / 9 If ^ 0 * 9 1? y U » / 

’ n , / 

c , ■■ . , I i, ( n i. /.(.,», i.i *'..5 f f ), r- 1 , ;n 

'! ■-' ': '1 t, D ;/.ti) ,0.1.1 ..3( rj , J = 1 ,J) 

''ll,!--' t(.ij, '!''7Uj,0,...nn,.l.s:l,.3) 

- .. r.s.',. »‘. ( { r- { / ) / K, ' , -ifv^A't.ih <:rfE;re') ■ 

I \ .,(/// ^ , '.''■■•'■•i '.t '-A iO), Cf-clKKV) 

/ , I . I' 

' ,„i !'!!.'-> A j (/,(, , 4 ) 

f’(<J i'l '»./■!, TJ 'if I) ,1 = 1,3) 

TvPr-/. > 3, 1 U ;.?/ 

i.-Tn,!.'!) ,Tf.*-) = l ■tf.dOU .,1) 

H'-} Cum I fit t 

DU M T=7,C,? 

c=r/2 

T-I'ii..;hv,1 ,3 )=(•'> 1KK(.1'-1. ,UMon,-!UM 
1 ."'’'uM I '+ 1 , H = ' I ( iil ^flC'i'CK f I , I ) 

,1 3 r.P’ C'l -lUti 

M I .U M =3 ,f> , '.i 

I ' = 1 ! -' / ' ■ 

1 ■;i; M C i l (X.) ♦tTwTfe.RlTl-i ,2)-!1r.'ftM2tL)tIN I’KaCll ,1 ) ) 

X .. Vi;,« a Ui , 2 )*‘Xu;‘ 1 ( !, ) 3=1 NTER ( 1 1 , 2 ) +!J!JP'2 CL) ♦INTERC 1 1 , 1 3 
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